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SUMMARY 

The stress and displacement distribution in a layered medium is 
found by means of transfer matrices. The surface loading exhibits axial 
symmetry, and each layer is of infinite extent in the horizontal direction, 
of constant depth, and is considered to be linearly elastic, homogeneous, 
and isotropic. The method developed has the built-in advantage of enforc- 
ing interface continuity conditions automatically. Its application to 
layered composites shows the flexibility with which it predicts the local 
as well as the global response of the medium. 

INTRODUCTION 

Recently, this writer developed a transfer matrix approach to various 
problems in mechanics by combining the method of initial functions due to 
Vlasov (ref. 1), with the integral transform method developed by Sneddon 
(ref. 2) . 

The method employed by this writer consists in applying the state 
space approach, which has been used extensively to analyze linear systems 
in various areas of systems engineering, such as modern control theory 
(ref. 3), to the field of elas t oinechanics . 

The topics so far analyzed through this approach cover two-dimensional 
elastostatics (ref. 4), one- dimensional elastodynamics (ref. 5), applica- 
tion to a typical elasticity problem (ref. 6), examination of the basic 
foundation of the theory (ref. 7), application to numerical integration of 
equations of motion to predict dynamic response (ref. 8), heat conduction 
(ref. 9), boundary value problems (ref. 10) and earthquake engineering with 
emphasis on soil-structure interaction (ref. 11). Additional references 
pertaining to each topic considered will be found in the references cited 
above and will not be repeated here. 
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This paper extends the work described in (ref. 4) which was restricted 
to a plane stress (or plane strain) problem, to a three-dimensional one with 
axially symmetric loading. The motivation for considering the present 
approach is to develop a flexible method for the analysis of layered media 
subjected for instance to concentrated loads, ranging from classical problems 
in soil mechanics, to the prediction of impulsive response of laminated com- 
posites. In the latter case inertial effects must be included. 

The main advantage of the method is due to the fact that continuity 
of stresses and displacements at interfaces is automatically satisfied. 
Therefore, upon determination of the missing initial displacements from 
boundary conditions, the field quantities can be determined upon multipli- 
cation of the initial state vector by the chain of layer transfer matrices 
by the field matrix of the layer of interest. A Hankel inversion gives the 
actual field quantities. 

In contrast, the classical formulation requires the construction of a 
transformed Airy stress function that contains four arbitrary parameters 
per layer, thus producing a total of 4n equations in 4n unknowns for a 
medium of n layers. These are determined by enforcing the continuity of 
stresses and displacements across each interface, which yields 4(n-l) con- 
ditions to which the four boundary conditions are added. 

DERIVATION OF THE TRANSFER MATRIX 


The equations governing the state of stress of an axially symmetric, 
homogeneous, isotropic, linearly elastic solid, are given by the equi- 
librium equations 
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The four stresses given by equations (2) are functions of the partial 
derivatives of two displacements only? it follows that two of these stresses 
can be eliminated. 


For reasons of convenience, Oq and a r are chosen for this purpose. 
Upon substitution of equations (2a) and (2b) into equation (la), the latter 
can be rewritten as 
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Elimination of the mixed derivative between equations (3) and (4) 
results in the relation 
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Consider a semi- infinite elastic medium which extends to infinity in 
the r-direction as shown in figure 1. The medium is loaded by an axially 
symmetric load as shown. Under the circumstances, taking Hankel Transforms 
of order one of equations (5) and (2d) , and of order zero of equations (2c) 
and (lb), results in the system of equations cast in matrix form as follows: 
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where the subscripts indicate the order of the Hankel transform. Equation 
(6) can be integrated by considering the column^vector of transformed 
stresses and displacements as the state vector X(£,z) , and rewriting it as 


{xa, Z )> = [A(g) ] {xa, Z )} (7) 

As shown in (ref. 4) s equation (7) can be integrated to yield 

{£a,z)} = exp[zA(?)] {X( 0 )} ( 8 ) 

where the matrix exponential has to be evaluated explicitly. The char- 
acteristic roots of the determinant associated with the matrix A(£) are the 
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double roots + £, identical to the result obtained in ref. 4 . Therefore, 
the results are analogous to those obtained in that paper, in which it is 
shown that 


where 
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Upon substitution of these values into equation (9), the transfer 
matrix is obtained, and equation ( 8 ) gives, in turn, the state vector which 
consists of the transformed stresses and displacements at an arbitrary depth 
in the field. The details pertaining to the evaluation of the transfer 
matrix are given in the Appendix. The results can be summarized in matrix 
form as 
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where the influence functions mapping the initial field quantities into 
those at an arbitrary depth in the field are given by 
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The actual physical quantities are then recovered through the inverse 
Hankel transform. 

APPLICATION TO LAYERED SYSTEM 


Consider a layered medium with perfect bonding along all interfaces 
as shown in figure 2. This implies the continuity of transformed stresses 
and displacements across each interface. In order to enforce this condi- 
tion, the first two entries of the state vector which appear in equation 
( 11 ) are divided by the shear modulus, to produce a new state vector con- 
sisting of transformed stresses and displacements. The elements of the new 
matrix G become Gp = Lp/pj G^ = Lp/p; G 2 3 = L 23 /P; G 24 = I '24^» 

G^ = pL^; G 32 = ,jL 32’ G 41 = ^41’ and G 42 = yL 42* 11116 remainin § elements 

of the G matrix are identical to the corresponding elements of the L matrix. 


The modified equation (11) can now be written in contracted form as 

'*?(£,*)}- [G(X,p,z,5)l (Y(5,0)} (13) 

Applying equation (13) to each interface in turn, in the sequence shown in 
figure 2, leads to 

(Y(?,h )} = [G(X ,p ,h ,C)--*G(A ,p ,h , 5 )] {Y(C,0)} (14) 

n n n n 111 

in which the missing initial conditions are determined from boundary con- 
ditions. Equation (14) then describes the overall response of the layered 
system. 

Local information consisting of state vectors at interfaces can now be 
obtained by terminating the matrix multiplication indicated by equation (14) 
at the appropriate interface. These relations are shown by the block 
diagrams shown in figures 3 and 4. 


The state vector in any arbitrary layer m can now be found by the 
relation 


m-1 

{Y(?,z)} = [G(A ,p ,z,Q II G(X ,p ,h.,C)]{Y(S,0)} (15) 

mm ±=1 x x 1 


in which the z coordinate is the local depth within the layer m, ranging 
from zero to h^. The actual stresses and displacements are given by the 
inverse Hankel transformation of the state vector. 


CONCLUDING REMARKS 

In this paper, a transfer matrix method to determine the response of 
a layered medium subjected to an axially symmetric loading has been 
presented. 


The matrix formulation shows that the need for matching interface 
conditions explicitly is avoided by imposing the continuity of the state 
vector across each interface. This is accomplished through the continued 
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multiplication of layer transfer matrices. Therefore, the size of the trans- 
fer matrix remains four by four, and is independent of the number of layers 
contained in the medium. This is the main conceptual as well as computa- 
tional advantage of the proposed method. 


APPENDIX 



The transfer matrix is given by the expression exp(zA) = a^I + a^A + 

3 2 3 

+ a^A , in which the matrices A,A , and A are given by 
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and the coefficients aQ,a^a 2 » and a^ are given by the set of 
relations (10) . The elements of the matrix exponential are given 
explicitly by the expressions (12) . 
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